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Various cooling methods, including a recently introduced one which smoothes out only quantum fluctuations larger than 
a given threshold, are applied to the study of topology in 2d CP'^"^ models. A critical comparison of their properties is 
performed. 



Two-dimensional CP^~^ models play an impor- 
tant role in quantum field theory because they share 
many properties with QCD. In particular, they pos- 
sess instanton classical solutions and one can define 
topological charge and susceptibility in analogy with 
QCD. CP'^^^ models represent therefore a useful 
theoretical laboratory to investigate numerical meth- 
ods to be eventually applied to study the topology in 
QCD. One of the most powerful tools for the study of 
the topological structure of the vacuum is the "cool- 
ing" method. It consists in measuring the topolog- 
ically relevant quantities on the ensemble of lattice 
configurations obtained by replacing each equilibrium 
configuration by the one resulting after a sequence of 
local minimizations of the action. 

The aim of this work is to get insight into a new 
cooling method first adopted in Ref. [IJ, by compar- 
ing it with the "standard" EJ and with its "controlled" 
version adopted by the Pisa group f3l. 

We chose the standard discretization for the action 
ofthe2dCP^-i model: 

— ~N [3 ^ ^ (2n+/j^nA„^^ + Z„Z„-|„^A„.p — 2) (1) 

where z„ is an iV— component complex scalar field 
with Zn- Zn = 1, A„^p is a U(l) gauge field satisfying 
An,;^A„.p ~ 1 and (3 ~ l/{Ng), with g the lattice 
coupling. We used the standard action both in Monte 
Carlo simulations and in the cooling instead of any 



improved lattice action, since we wanted to test the 
cooling techniques in a situation where cutoff effects 
are large. The lattice topological charge density was 
defined as 

9^ (a;) = 5] e^.Tr [P{x)A,P{x)A,Pix)] , (2) 

fiiy 

with 

Pij{x) = Zi{x)zj{x), (3) 
giving for the lattice topological susceptibility 

X^^(^g^(.)g^(0)\ = -^((Q^)') , (4) 

where = J^x^'^i^) '^^'^ ^ lattice size. 

The cooling algorithm consists in assigning to each 
lattice variable z„ and A„.p a new value 2;"°* and 
A^™ which locally minimizes the action, keeping all 
other variables fixed. In the "standard cooling" these 
replacements are unconstrained. We will call "new 
cooling" the one for which the replacements are done 
only if the angle a between the new and the old field 
variables, is larger than a given value 6, and "Pisa 
cooling" the one for which the local minimization is 
performed with the constraint a < Spisa,. Notice that 
between the "Pisa cooling" and the "new cooling" 
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Figure 1. vs p/a under the three different cool- 
ing procedures. The instanton size is determined by 
^Lmax ^ 2/{p/af, (a is the lattice spacing).The 
three curves faU exactly on top of each other 

there is a substantial difference: while "Pisa cooling" 
acts first on the smoother fluctuations, the "new cool- 
ing" performs local minimizations only if these fluc- 
tuations are larger than a given threshold. Moreover 
the "new cooling" automatically stops when there are 
no more fluctuations beyond the threshold. It should 
be pointed out that any cooling procedure causes a 
partial loss of the topological content of the cooled 
configuration (namely "small instantons"). However 
this loss occurs at a fixed scale in lattice units and thus 
vanishes in the continuum limit, unless the instantons 
distribution is ultraviolet singular (as in the case of 
CPi). 

We considered first an "artificial" 1-instanton con- 
figuration discretized on the lattice. Although the 
three different cooling procedures act differently, 
since they start to deform the configuration in dif- 
ferent regions, the curves for and under the 
three coolings as a function of the instanton size p fall 
on top of each other (see Fig.Q. We determined the 
topological susceptibility using the "new cooling" and 
compared the results with those from the "field theo- 
retical method" [4J. In the field theoretical approach 
one has 

X^(/3) = a'Z{(3fx^°^' + M{(3) + 0{a^) 
M{P) - a^A{P){S{x))np+PW{I), (5) 



where Z{P) is a finite multiplicative renormaUza- 
tion of the discretized topological charge density |0, 
while M{(3) is an additive renormalization contain- 
ing mixing to operators of equal or lower dimension 
and same quantum numbers, namely the action den- 
sity and the identity operator. On the lattice 
is measured during the Monte-Carlo simulation, and 
^cont extracted by subtracting the renormalizations, 
which are computed non-perturbatively on the lattice 
by means of the "heating method" f6l (Z and P can 
be also computed perturbatively [7J). The "field theo- 
retical method" can be improved by using a smeared 
topological charge density operator fsl, built from the 
standard operator by replacing the fields z and 
with smeared fields 2;^™'=^'' and A™^''' . In this way 
the renormalizations are strongly reduced and a much 
better accuracy is achieved for ^'^Q"*. 

We performed numerical simulations for = 4, 
10, 21. The simulation algorithm is a mixture of 4 
microcanonical updates and 1 over-heat bath. For 
each simulation we collected lOOK configurations af- 
ter lOK thermalization updating steps. We used sev- 
eral values of the 5 parameter in the "new cooling" 
algorithm, while for the "Pisa cooling" we considered 
5pis& = 0.2. To set the scale we have taken the cor- 
relation length defined as the second moment of 
the correlation function (YrP{x) P{0)). In Fig.|2lwe 
show the behavior of on configurations cooled by 
the "new cooling" with different (5's, while in Fig. [S] 
we plot x'^°'^^(,G f^'^ values of S which correspond 
to the peak region in x^ in Fig- IZl ™d compare the 
results with the field theoretical determination (with 
and without smearing). There is consistency between 
the two determinations for all the considered values of 
6, which correspond to those for which the smoothing 
of the configuration is better. We have also caUbrated 
the three cooling techniques (number of cooling steps 
for the standard and Pisa cooling versus 5 for the new 
cooling) in order that the average energy on the en- 
semble cooled in the three different ways is the same. 
Then, comparing by eye several thermal configura- 
tions obtained after equivalent amounts of the three 
coolings, we have observed that the distributions and 
the shape of the instanton bumps is roughly the same. 
Also the values of x'^°"'?g' obtained after "equiva- 
lent" coolings are in agreement within the statistical 
errors. The only measurement where we could see a 
discrepancy is that of the "shell" correlation function 
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Figure 2. as a function of 6 at (3 — 0.70 on a 24^ 
lattice for CP9. 

of the topological charge density ((Zl('')'?l(0)) (see 
Fig. 13, which is slightly larger in the case of the new 
cooling method in the short distance region, with re- 
spect to the standard and Pisa coolings. 

Our conclusion is that, except for the last observa- 
tion which deserves further study, there is no appre- 
ciable difference between the three types of cooling 
we have investigated. 
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Figure 3. x'^^^Cg forCP^ for several S values (L/^ : 
10 for all the points). 
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Figure 4. "Shell" correlation function of the lattice 
topological charge density (r is the distance between 
two lattice sites). 



